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I have measured the early stages of the growth of branched 
metal aggregates formed by electrochemical deposition in very 
thin layers. The growth rate of spatial Fourier modes is de- 
scribed qualitatively by the results of a linear stability anal- 
ysis [D.P. Barkey, R.H. Muller, and C.W. Tobias, J. Elec- 
trochem. Soc. 136, 2207 (1989)]. The maximum growth 
rate is proportional to (I/c) s where I is the current through 
the electrochemical cell, c the electrolyte concentration, and 
S = 1.37±0.08. Differences between my results and the theo- 
retical predictions suggest that electroconvection in the elec- 
trolyte has a large influence on the instability leading to ram- 
ified growth. 
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Rough surfaces develop in many systems in which the 
surface grows under far-from-equilibrium conditions , 
due to the fact that the smooth surface is unstable to 
the growth of perturbations over a range of wave num- 
bers. This morphological instability was first treated 
by Mullins and Sekerka (MS) in the context of solidi- 
fying alloys H . Metal aggregates grown by electrochem- 
ical deposition (ECD) in quasi-two-dimensional geome- 
tries display a variety of qualitatively different branched 
growth morphologies and there has been consider- 

able study of the basic instability of the straight elec- 
trode which leads to the initial development of branched 
growth §-|l|. 

To understand this instability qualitatively, assume 
that the cations are transported by migration in an ap- 
plied electric field to the cathode, where they are de- 
posited. The electric field will be uniform along the 
length of a straight cathode, but will be enhanced near 
the tip of a small bump. This will lead to an increased 
current of cations, and so a higher deposition rate, at 
the tip of the bump. Thus the bump will grow expo- 
nentially. A similar argument applies if the transport is 
diffusive. At very high spatial frequencies, however, the 
straight interface is stabilized by surface tension. This 
picture leads to an instability of the straight electrode 
analogous to the Mullins-Sekerka instability — pertur- 
bations within a band of wave numbers k from k = up 
to a cutoff wave number k c will grow, while perturbations 
with k > k c will be damped out. The exponential growth 
rate (3 of spatial modes with wave number k has the form 



14fl studied growth by 
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with k c = r^ 1 / 2 . 

Experimentally, Kahanda et al. 
electrochemical deposition in the 

rents. In this case Jl4j], the deposition of metal onto the 
cathode is limited by the activation of cations in the dou- 
ble layer. Kahanda et al. used a Fourier analysis tech- 
nique to study the growing aggregate front. They deter- 
mined the growth rate of the Fourier modes as a function 
of wave number and found semi-quantitative agreement 
with the predictions of MS theory. 

At higher currents, however, the deposition rate is gov- 
erned by the transport of ions in the electrolyte. A 
number of electrochemical processes contribute to the 
ion transport and to the current distribution near the 
cathode, resulting in modifications to Eq. |l| ]7| p|, pd| , ^3| . 
Barkey, Muller, and Tobias (BMT) jyj have performed 
a linear stability analysis of a planar electrode in three- 
dimensional ECD, and, at least conceptually, the results 
of their work should carry over to the two-dimensional 
case studied here. They assume mass transport to be 
due solely to diffusion. Surface tension manifests itself 
as a reduction in the binding energy of deposited atoms 
on a curved surface relative to a flat surface (referred to 
in BMT as the capillary potential shift). The effects of 
electrode reactions (kinetic overpotential) and of charge 
separation in the diffusion layer close to the cathode (con- 
centration overpotential) both act to smooth out varia- 
tions in the current density, and thus reduce the growth 
rate of perturbations. BMT found (see Eq. (41) of Ref. 
|nj) the growth rate (3{k) to have the form 



qk(l — rk 2 ) 
(1 + sk) ' 
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The coefficients q, r, and s depend on the properties of 
the electrolyte and the deposited metal, the ion concen- 
tration, and the current Jl5[ . Roughly speaking, q in- 
corporates the destabilizing effects of mass transport, qr 
involves the stabilizing effects of surface tension, and s is 
a combination of terms involving the kinetic overpoten- 
tial and the concentration overpotential. As above, there 
is a range of unstable wave numbers, < k < r -1 / 2 ; the 
MS result is recovered in the limit that s — > 0. 

In this Communication, I present results from experi- 
ments on the very early stages of the ECD of aggregates 
of metallic copper from solutions of CuSO.4, in the regime 
where the growth is limited by mass transport in the 
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electrolyte. Usin g an alysis techniques similar to those 
of Kahanda et al. B], I investigate the growth rate (3{k) 
and find it to be qualitatively described by the dispersion 
relation of BMT (Eq. §) and not by Eq. 0. 

In my experiments [lq], two copper foil electrodes 5.1 
cm long by 0.025 cm thick, separated by roughly 23 mm, 
were sandwiched between two 5.1 cm square by 0.6 cm 
thick glass plates and clamped together. The space be- 
tween the electrodes was filled with aqueous solutions of 
CuSC>4 with concentrations c in the range 0.02 M < c < 
0.5 M. A constant current I of from 0.06 mA to 10 mA 
was passed through the cell, and a branched copper ag- 
gregate formed at the cathode. The cell was illuminated 
from below by a diffuse white light source, and imaged 
from above with a ccd video microscope. Individual im- 
ages of the aggregate were captured and digitized by a 
frame grabber in a personal computer, and the growth 
was also recorded on video tape. The spatial resolution 
of the digitized images was typically 23 /im/pixel. 

My data analysis was similar to that of Kahanda et 
al. JL4J The edge of the aggregate (the interface) was lo- 
cated by thresholding the digitized image. Pixels with 
an intensity lower than a chosen value were deemed to 
be on the aggregate. Since the aggregate is rough, the 
interface determined in this way, p(x) 7 will in general 
not be a single- valued function of the coordinate x along 
the length of the cathode. For my analysis, I formed a 
single- valued approximation to the interface, pi(x), by 
taking the point on the interface furthest from the orig- 
inal cathode position, i.e., the highest point, at each x 
value. For a very heavily branched aggregate with many 
"overhangs," pi(x) is not a good representation of the 
true interface. However, for the early times of interest 
here, the aggregate, although rough, has few overhangs, 
and pi(x) is very similar to p(x). The single- valued in- 
terface function p\ (x) is then Fourier transformed to g ive 
a spatial Fourier power spectrum of the interface |17f| . 

Fig. ^ shows the interface functions determined at two 
different times for a particular run. Curve (a) of Fig. 
[j] shows p(x) determined from a digitized video image 
recorded 225 s after the start of the run. The aggre- 
gate is rough, but, at least within the resolution of my 
imaging system, it does not yet contain a large number 
of overhangs. Curve (b) of Fig. |l| is the single-valued 
function p\{x) extracted from curve (a); it is very sim- 
ilar to p(x). Curve (c) shows p(x) at a later time, 385 
s after the start of the run. By this time the aggre- 
gate is more strongly branched and displays substantially 
more overhangs. Curve (d) is the corresponding single- 
valued function p\(x). In this case there are significant 
differences between p\(x) and p{x) and Fourier analysis 
of p\ (x) does not provide meaningful information about 
the true interface. 

By Fourier analysing the single-valued interface func- 
tions obtained from a time-sequence of images, one can 
extract the Fourier power |v4(/c)| 2 of a mode of wave num- 
ber A: as a function of time. Fig. |] shows |^4(fc)| 2 vs. time 
for three different values of k, from the run illustrated in 



Fig. [|. In this particular run the growth at the cath- 
ode first became visible about 100 s after the start of the 
run. The Fourier power grows exponentially with time 
between this time, and about 230 s. The dotted lines in 
Fig. U are fits of the data to a growing exponential over 
this time range. The growth rate depends on k. At times 
beyond the region of exponential growth, the growth rate 
generally decreases. In this regime, however, p\{x) is not 
a good approximation to the true interface and quanti- 
tative analysis is not possible using this technique. The 
interface shown in curves (a) and (b) of Fig. [l] corre- 
sponds to a time near the end of the exponential growth 
phase of Fig. |2|, while that of curves (c) and (d) of Fig. 
[I] is well beyond the end of this phase. 

Fig. H shows the exponential growth rate (3 of the 
Fourier amplitude \A(k)\ as a function of k. The dotted 
curve is a least-squares fit to Eq. [j], which does not de- 
scribe the data well. Rather, the growth rate increases 
more rapidly for small k, displays a weak maximum, and 
then decreases more slowly for larger k. The solid curve 
plotted in Fig. || is a least-squares fit of the data to a dis- 
persion relation with the form of Eq. ||. The fit is quite 
good, although, as in Fig. |L most data sets show per- 
haps a bit more of a peak in the data than the theoretical 
curve would indicate. Although fits of the experimental 
data to Eq. || looked satisfactory to the eye, the uncer- 
tainties in the parameters were in many cases quite large. 
This was due to the fairly large scatter in the (3{k) data; 
that in the data shown in Fig. |^ is typical. 

From fits of Eq. || to the data from a number of runs, 
we determined the growth rate fimax of the fastest grow- 
ing mode. Fig. |J(a) is a plot of (3 m ax against the cell 
current / for two values of the concentration c. Although 
the scatter in the data is considerably larger than the er- 
ror bars, the trend of the data is well described by the 
power law (3 m ax oc I Sl with an exponent Si = 1.52 ± 0.19 
for c = 0.02 M and S t = 1.27 ± 0.06 for c = 0.1 M. Fig. 
^(b) shows Pmax as a function of c for a fixed current. In 
this case a power law fit to the data gives flmax °c c~ d2 , 
with §2 = 1-40 ± 0.16. The similarity between the values 
of the exponents Si and 5% suggests that Pmax may be 
a function of the ratio I/c alone. In Fig. ||(c) the same 
three data sets are plotted as a function of I/c. With 
some scatter, the data collapse onto a single power-law 
function given by (3 ma x = (5.2± 1.2) x lO-^I/c) 1 - 37 ^ - 08 
s~l. 

To compare my results with the theoretical predictions 
|]llf , I make the assumption that my measurements were 
performed at currents such that i < i^, where i is the 
current density and %l the limiting current density p8| , 
and use bulk values for the electrolyte properties involved 
in the expressions for the coefficients q, r, and s |Tq| . The 
theory then predicts that the coefficients will depend on 
current and concentration as q oc /, r tx c/J, and s oc c/I. 
On the other hand, the values of q obtained from fits to 
the experimental data tend to increase with I at constant 
c, but decrease with c at constant /; the fitted values of 
r show no systematic variation with either / or c and in 
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fact are constant to within a factor of two over the range 
of conditions studied; and the values of s, while they have 
large error bars, also show no systematic variation with 
either I or c. Numerically, for c = 0.1 M and I = 2mA, 
the value of q obtained from fits to the experimental data 
is three orders of magnitude larger than predicted, r is 
a factor of 20 larger than predicted, while s is approxi- 
amtely equal to the predicted value. Under the same as- 
sumptions, the theory predicts that Pmax °c I 2 /c, while 
Fig. H(c) indicates that, experimentally, fimax behaves 
as a power law in 7/c. 

The fact that the shape of the experimental dispersion 
relation can be fitted by Eq. suggests that the basic 
ingredients of the theory — that is, that the instability 
is driven by mass transport, stabilized by surface ten- 
sion, and modified by electrochemical processes — are 
correct. However, the quantitative differences between 
experiment and theory outlined above suggest that ei- 
ther the assumptions made by me in calculating the co- 
efficients are invalid, or that effects not accounted for in 
the theory play an important role in these experiments. 

Both q and s are predicted to diverge when i — [ p"5[ . 
In calculating the coefficients, I assumed that i <C il- A 
rough estimate of the limiting current for my cells JlS| ] 
gives II ~ 5 mA for c = 0.1 M, similar to the highest 
currents used at that concentration. However, to get a 
thousand-fold increase in q, one would need to be within 
0.1% of iz- I find larger than expected values of q over a 
large range of currents, and the assumption that i <C ih 
is safe over most of that range. On the other hand, 
it is quite likely that the use of the bulk value of the 
electrolyte conductivity in evaluating the coefficients is a 
poor approximation, since in an unsupported electrolyte 
a diffusion layer depleted of ions, and with a correspond- 
ingly low conductivity, develops near the cathode. How- 
ever, using a conductivity of zero would decrease q by 
roughly a factor of two, while r and s would go to zero 
p5[ . Thus these assumptions cannot account for the ob- 
served discrepencies. 

The theory of BMT assumes that ion transport is due 
to diffusion only. In an unsupported binary electrolyte, 
migration will also be important. This can be accounted 
for by a correction factor which would result in a decrease 
of the theoretical value of q by about 20% Q . 

The most likely explanation for the differences be- 
tween my experimental results and the theoretical pre- 
dictions is the fact that high electric fields in the dif- 
fusion layer cause electroconvectivc flow in quasi-two- 
dimensional ECD experiments |2(||2l|. This flow has an 
influence on ion transport and on the deposition process, 
and is not accounted for in the theory of BMT. Natural 
convection, driven by density gradients in the electrolyte 
[p3|-p6[ may similarly play a role. 

In summary, I have measured the growth rate of spatial 
modes as a function of wave number during the early 
stages of the growth of metal aggregates in thin-layer 
electrochemical deposition. The dispersion relation is not 
well described by a MS-type theory m , but can be fitted 



to the form predicted by the stability analysis of Ref. [[Ll| . 
Significant quantitative differences between my results 
and the predictions of the theory of Ref. |ll[] suggest that 
convection in the electrolyte driven by electric fields or by 
density gradients modifies the ion transport substantially 
and plays an important role in determining the growth 
rate. A more detailed account of this work will appear 
elsewhere p9| . 

This research was supported by the Natural Sciences 
and Engineering Research Council of Canada. I am 
grateful to D. Barkey, G. Marshall, G. White and S. Mor- 
ris for helpful discussions. 
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FIG. 1. The edge of the metal aggregate growing on the 
cathode for a run with c = 0.1 M and I = 1.4225 mA. (a) 
The interface function p(x) at time t = 225 s, near the end 
of the exponential growth regime seen in Fig. | below, (b) 
The single- valued function pi(x) corresponding to curve (a), 
(c) p(x) at t = 385 s, well beyond the end of the exponential 
regime, (d) The single- valued function pi(x) corresponding 
to curve (c). The solid bars have length 1mm. 

FIG. 2. The Fourier power as a function of time for three 
values of wave number k, for the same run as shown in Fig. 
|IJ Triangles: k — 27.2 mm -1 , circles: k = 75.3 mm -1 , dia- 
monds: k = 113 mm -1 . The lines are fits to the data in the 
exponential growth regime. 

FIG. 3. The growth rate f3(k) of the Fourier amplitude 
\A\, for the same run as the previous figures. The dotted 
curve is a fit to Eq. [j], and does not describe the data well. 
The solid curve is a fit to the form predicted by BMT, Eq. ^[ 

FIG. 4. (a) The maximum growth rate fimax as a function 
of the current I for fixed concentration c. Squares: c = 0.02 
M; triangles: c = 0.10 M. The lines are fits of the data to 
power laws in /, as discussed in the text, (b) Pmax as a 
function of c for I = 2.00 mA. The dashed line is a fit to a 
power law. (c) The same data as in (a) and (b) plotted as a 
function of 7/c. The symbols are as in (a) and (b), and the 
line is a fit of all the data to a power law. 
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